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The LDA+DMFT approach merges conventional band structure theory in the local den- 
sity approximation (LDA) with a state-of-the-art many-body technique, the dynamical 
mean-field theory (DMFT). This new computational scheme has recently become a pow- 
erful tool for ab initio investigations of real materials with strong electronic correlations. 
In this paper an introduction to the basic ideas and the set-up of the LDA+DMFT 
approach is given. Results for the photoemission spectra of the transition metal ox- 
ide Lai—xSr^TiOa, obtained by solving the DMFT-equations by quantum Monte-Carlo 
(QMC) simulations, are presented and are found to be in very good agreement with 
experiment. The numerically exact DMFT(QMC) solution is compared with results ob- 
tained by two approximative solutions, i.e., the iterative perturbation theory and the 
non-crossing approximation. 
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1. Introduction 

It is an established fact of solid state physics that many electronic properties of 
matter are well described by the purely electronic Hamiltonian 
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where the crystal lattice enters only through an ionic potential. The applicability 
of this approach may be justified by the validity of the Born and Oppenheimer 
approximations. Here, $ + (r,a) and ^(r, a) are field operators that create and 
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annihilate an electron at position r with spin a, A is the Laplace operator, m e the 
electron mass, e the electron charge, and 

Mo„(r) = - e 2 £-^- and Ke(r-r') = £ £ * (2) 

denote the one-particle potential due to all ions z with charge at given positions 
Ri, and the electron-electron interaction, respectively. 

While the ab initio Hamiltonian ([!]) is easy to write down it is impossible to 
solve, even numerically, if more than a few electrons are involved. The reason 
for this is the electron-electron interaction which correlates every electron with all 
others. Therefore, one either needs to make substantial approximations to deal with 
the Hamiltonian ([[]), or replace it by a strongly simplified model Hamiltonian. At 
present these different strategies for the investigation of the electronic properties 
of solids are applied by two largely separate groups: the density functional theory 
(DFT) and the many-body community. It is known for a long time already that 
DFT, together with its local density approximation (LDA), is a highly successful 
technique for the calculation of the electronic structure of many real materialscl. 
However, for strongly correlated materials, i.e., d- and /-electron system which 
have a Coulomb interaction comparable to the band- width, DFT/LDA is seriously 
restricted in its accuracy and reliability. Here, the model Hamiltonian approach 
is more general and powerful since there exist systematic theoretical techniques 
to investigate the many-electron problem with increasing accuracy. These many- 
body techniques allow to describe qualitative tendencies and understand the basic 
mechanism of various physical phenomena. At the same time the model Hamiltonian 
approach is seriously restricted in its ability to make quantitative predictions since 
the input parameters are not accurately known and hence need to be adjusted. One 
of the most successful techniques in this respect is the dynamical mean-field theory 
(DMFT) - a non-perturbative approach to ateaagly correlated electron systems 
which was developed during the past decadeHu'Ercl The LDA+DMFT approach, 
which was first formulated by Anisimov et al.0, combines the strength of DFT/LDA 
to describe the weakly correlated part of the ab initio Hamiltonian (|l|), i.e., electrons 
in s- and p-orbitals as well as the long-range interaction of the d- and /-electrons, 
with the power of DMFT to describe the strong correlations induced by the local 
Coulomb interaction of the d- or /-electrons. 

Starting from the ab initio Hamiltonian (|l|), the LDA+DMFT approach is 
derived in section 2. As a particular example, the LDA+DMFT calculation for 
Lai_ x Sr x Ti03 is discussed in section 3. Furthermore, the LDA+DMFT spectrum 
is calculated by means of numerically exact quantum Monte Carlo (QMC) simu- 
lations, and is compared to that obtained within two approximations commonly 
employed to solve the DMFT, i.e., the iterated perturbation theory (IPT) and the 
non-crossing approximation (NC A) . A discussion of the LDA+DMFT approach and 
its future prospects (section 4) closes the presentation. 
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2. The LDA+DMFT approach 



2.1. Density functional theory 

The fundamental theorem of DFT by Hohenberg and KohiJ] states that the ground 
state energy is a functional of the electron density which assumes its minimum at 
the ground state electron density. Following Levy!, this theorem is easily proved 
and the functional even constructed by taking the minimum (infimum) of the energy 
expectation value w.r.t. all (many-body) wave functions if(r±ai, ■ ■■'F'no'n) at a given 
electron number N which yield the electron density p(r): 

N 

E\p]=ini{{<p\H\<p) | (tp\ 1 £S(r-r i )\<p}=p(r)}. (3) 

t=i 

However, this construction is of no practical value since it actually requires the 
evaluation of the Hamiltonian (|l|). Only certain contributions like the Hartree 
energy ^Hartree [p] = 5/ d 3 r' d 3 r V cc (r — r') p(r')p(r) and the energy of the ionic 
potential Ei on [p] = J d 3 r Vi on (r) p(r) can be expressed directly in terms of the 
electron density. This leads to 

E [p] = Skin [p] + Son [p] + S H artrcc [p] + S xc [p] , (4) 

where Sk; n [p] denotes the kinetic energy, and E xc [p] is the unknown exchange and 
correlation term which contains the energy of the electron-electron interaction ex- 
cept for the Hartree term. Hence all the difficulties of the many-body problem have 
been transfered into E xc [p\. While the kinetic energy Skin cannot be expressed 
explicitly in terms of the electron density one can employ a trick to determine it. 
Instead of minimizing E[p] w.r.t. p one minimizes it w.r.t. a set of one-particle 
wave functions ifi related to p via 

N 

p(r)=5>,(r)| 2 . (5) 

»=1 

To guarantee the normalization of ipi, the Lagrange parameters £i are introduced 
such that the variation 5{S[p] + e s [l — J d 3 r|( ( Cj(r)| 2 ]}/(5i J 9j(r) = yields the Kohn- 
Sham0 equations: 



l-A + V ion (r)+ J d 3 r'p(r')V ee (r-r') + 



2ni, J ' ' Sp(r) 

These equations have the same form as a one-particle Schrodinger equation which, 
a posteriori, justifies to calculate the kinetic energy by means of the one-particle 
wave-function ansatz. The kinetic energy at the ground state density is, then, given 
by i?kin[/3min] = — YliLi (fi \fo 2 A/ (2m e ) \<pi) if the ipi are the self-consistent (spin- 
degenerate) solutions of Eqs. (||) and (||) with lowest "energy" e,. Note, however, 
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that the one-particle potential of Eq. (^), i.e., 

Mo„(r) + Jd 3 rV(r')y cc (r-r') + (7) 

is only an auxiliary potential which artificially arises in the approach to minimize 
E[p\. Thus, the wave functions <pi and the Lagrange parameters Si have no physical 
meaning at this point. 



2.2. Local density approximation 

So far no approximations have been employed since the difficulty of the many-body 
problem was only transferred to the unknown functional E xc [p\. For this term the 
local-density approximation (LDA) which approximates the functional E xc [p] by a 
function that depends on the local density only, i.e., 

E*cM- / d 3 r E^ A (p(r)), (8) 

was found to be unexpectedly successful. Here, E x J ^ A {p{v)) is usually calculated 
from the Hartree-Fock solution or the numerical simulation of the jellium problem 
which is defined by V\ on {r) = constll3. 

In principle DFT/LDA only allows one to calculate static properties like the 
ground state energy or its derivatives. However, one of the major applications of 
LDA is the calculation of band structures. To this end, the Lagrange parameters e, 
are interpreted as the physical (one-particle) energies of the system under consider- 
ation. Since the true ground-state is not a simple one-particle wave-function, this 
is a further approximation beyond DFT. Actually, this approximation corresponds 
to the replacement of the Hamiltonian ([!]) by 

H LDA = fd 3 r*+(r,a) -A-& + V ion {v) + j d 3 r' p^V^r-v') 

U L 

, SE^[p] 



Sp{r) 



*(r,£r). (9) 



For practical calculations one needs to expand the field operators w.r.t. a basis 
<&Um, e.g., a linearized muffin-tin orbital (LMTO)0 basis (here i denotes lattice 
sites; I and m are orbital indices). In this basis, 



*+(r,a) = E^l-Wr) (10) 

ilm 

such that the Hamiltonian (^) reads 

^LDA = (Silrnjl'rn' £ilm fi>Um + hlm,jl' m' &U m &jV to') * (H) 

ilm, j I'm* ,<j 
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Here, nf lm — c ilm cf lm , 



i f z A r 5E \o\ \ 

tumjl'm' = (*«m ~^+ V ™( r ) + / d rV(r / )y c(r-r / ) + gfr) &j Vmf ) (12) 



for iim ^ jl'm! and zero otherwise; e^m denotes the corresponding diagonal part. 

As for static properties LDA band structure calculations are also often highly 
successful - but only for weakly correlated materialsa. Indeed, the self-consistent 
solution of the one-particle Hamiltonian -£Tlda ([ll]) together with Eq. (||) treats 
electronic correlations only rudimentarily. Consequently, LDA is not reliable when 
applied to correlated materials, and can even be completely wrong. For example, 
it predicts the antiferromagnetic insulator La2Cu04 to be a non-magnetic metaltj. 

2.3. Supplementing LDA with local Coulomb correlations 

Of prime importance for correlated materials are the local Coulomb interactions 
between d- and /-electrons on the same lattice site since these contributions are 
largest. This is due toJJie extensive overlap of these localized orbitals which results 
in strong correlationsEJ. To correct for these contributions, one can supplement 
the LDA Hamiltonian ([ll]) with the local Coulomb interaction XJ™ m , between the 
localized electrons (for which we assume i = id and I = Id)' 

-ffLDA+corrcl = -&LDA + - ^ U™ m ihii ma hii m i a > — #LDA ■ (13) 

i— id ,1= Id^raara' a 

Here, the prime on the sum indicates that at least two of the indices of an operator 
have to be different, and a term -Hlda i s substracted to avoid double-counting of 
those contributions of the local Coulomb interaction already contained in i?LDA- 
Since there does not exist a direct microscopic or diagrammatic link between the 
model Hamiltonian approach and LDA it is not possible to express H^da rigorously 
in terms of U and p. Guided by the observation that the LDA calculates the total 
energy of isolated atoms rather well, it was arguecO that the average energy .Elda 
corresponding to H^a * s we ^ approximated by the energy of the interaction term in 
the atomic limit. Hence, in the case of an orbital- and spin-independent U™ m , = U 
one may write 

^lda = \un d {n d - 1). (14) 

(For the corresponding equation including Hund's rule coupling see Rcf. I). Here, 
rid = J2 m n Hdm = ^2 m ( n n=idm) is the total number of interacting electrons. Since 
the one-electron LDA energies can be obtained from the derivatives of the total 
energy w.r.t. the occupation numbers of the corresponding states^ the one-electron 
energy level for the non-interacting states of (O) is obtained asO 



4 dm ■= j ~~ (Elda - ££ DA ) = e ildTn - U(n d - b (15) 
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where £n dm is defined in (|llj) and -Elda is the total energy calculated from i?LDA 
©■ 

This leads to a new Hamiltonian describing the non-interacting system 

-^LDA = ^ (^ilm,jl'm' £ ilm n ilm + tilm,jl'm'<%l m Cjl>m')i U^) 

ilm,jl'm f ,(T 

where e\ ldm is given by §U^) for the interacting orbitals and e\ lrn — eu m for the 
non-interacting orbitals. While it is not clear at present how to systematically 
subtract H^ DA one should note that the subtraction of a Hartree-type energy does 
not substantially affect the overall behavior of a strongly correlated paramagnetic 
metal in the vicinity of a Mott-Hubbard metal-insulator transition (see also Sec. 
2.6). 

In the following it is convenient to work in reciprocal space where the matrix 
elements of H®da are given by 

(-^LDA(k))gim,,j'2'm' = (HhT)AO<))qlm,q'l'm' $qlm,q' I' m' 3ql,q d l d U (rid ~ — ) . (17) 

Here, q is an index of the atom in the elementary unit cell, (-ffLDA(k)) 9 i m , g 'i' m ' is 
the matrix element of ([ll]) in k-space, and qd denotes the atoms with interacting 
orbitals in the unit cell. The non-interacting part, -ff^^A' supplemented with the 
local Coulomb interaction forms the (approximated) ab initio Hamiltonian for a 
particular material under investigation: 

-ffLDA+corrcl = -^LDA + J]] ^mm' Kilmer nilm' a' ■ (18) 

i—i d ,l—l d: mam' a 

To make use of this ab initio Hamiltonian it is still necessary to determine the 
Coulomb interaction U. To this end, one can calculate the LDA ground state 
energy for different numbers of interacting elecrons nd ("constrained LDA"EZl) and 
employ Eq. (|l4| ) whose second derivative w.r.t. nd yields U. However, one should 
keep in mind that, while the total LDA spectrum is rather insensitive to the choice 
of the basis, the calculation of U strongly depends on the shape of the orbitals 
which are considered to be interacting. Thus, an appropriate basis like LMTO is 
mandatory and, even so, some uncertainty in U remains. 

2.4. Dynamical mean-field theory 

The many-body extension of LDA, Eq. ([i"8|), was proposed by Anisimov et al0 
in the context of their LDA+U approach. Within LDA+U the Coulomb interac- 



tions of (18) are treated within the Hartree-Fock approximation; hence it does not 
contain true many-body physics. While LDA+U is successful in describing long- 
range ordered, insulating states of correlated electronic systems it fails to describe 
strongly correlated paramagnetic states. To go beyond LDA+U and capture the 
many-body nature of the electron-electron interaction, i.e., the frequency depen- 
dence of the self-energy, various approximation schemes have been proposed and 
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applied recentl^BBBEIl. One of the most promisingannroaches, first imple- 
mented by Anisimov et al.O, is to solve @ within DMFTQ'lii ("LDA+DMFT"). 
Of all extensions of LDA only the LDA+DMFT approach is presently able to de- 
scribe the physics of strongly correlated, paramagnetic metals with well-developed 
upper and lower Hubbard bands and a narrow quasiparticle peak, which is deter- 
mined by the vicinity of a Mott-Hubbard metal-insulator transition. 

During the last ten years DMFT has proved to be a successful apjproach to in- 
vestigate strongly correlated systems with local Coulomb interactionsQ. It becomes 
exact in the limit of high lattice coordination numbers! and preserves the dynamics 
of local interactions. Hence, it represents a dynamical mean-field approximation. 
In this non-perturbative approach the lattice problem is mapped onto an effective 
single-site problem which has to be determined self-consistently together with the 
k-integrated Dyson equation connecting the self energy E and the Green function 
G at frequency u>: 

Gqlm,q'l'm'(u) =— J d 3 k [u! Sqlm,q'l'm' ~ C^LDAO*))gJm.,g'i'm' 

_ Sql,q d l d ^qlm,q'Vm'{^)]~ ■ (19) 

Here, implies the inversion of the matrix with elements n (=qlm), n'(=q'l'm'), 

and integration extends over the Brillouin zone with volume Vb ■ 

The DMFT single-site problem depends on Q^ 1 = G^ 1 + E and is equivalent 
to an Anderson impurity problcmiH if its hybridization A(a->) satisfies Q (w) = 
ijj — J du>' A(ui') / (ui — uj'). The local one-particle Green function at a Matsubara 
frequency w„ = (2u + 1)tt//3, orbital index m (I = Id, q = qd), and spin a is given 
by the following functional integral over Grassmann variables ip and tp*: 

G° vm = -\\ V^\VWW vm r v * m ^ r ^ 1] - (20) 



Here, the single-site action A has the form 

E 



~\ E' / ^CWfi (r). (21) 

Due to its equivalence to an Anderson impurity problem a variety of approximative 



techniques have been emp 
perturbation theory (IPT) 



25j2a, m 



d to solve the DMFT equations, such as the i^t 
and the non-crossing approximation (NCA) 
as well as numerical techniques like quantum Monte-Carlo simulations (QMC)E 
exact diagonalization (ED)cJ'El, or numerical renormalization group (NRG)EJ. IPT 
is non-self-consistent second-order perturbation theory in U for the Anderson impu- 
rity problem (^p| ) at half-filling. It represents an ansatz that also yields the correct 
perturbational C/ 2 -term and the correct atomic limit for the self-energy off half- 
fillingpy. NCA is a resolvent perturbation theory in the hybridization parameter 
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A(cj) of the Anderson impurity problem. Thus, it is reliable if the Coulomb in- 
teraction U is large compared to the band-width. In essence, the QMC technique 



maps the interacting electron problem (20) onto a sum of non- interacting prob- 



lems by means of Hubbard -Straton ovi ch transformations and evaluates this sum by 
Monte-Carlo samplingtS. ED directly diagonalizes the Anderson impurity problem 
at a limited number of lattice sites. NRG first replaces the conduction band by 
a discrete set of states at DAT 11 (D: bandwidth; n = 0, ...,Af) and then diago- 
nalizes this problem iteratively with increasing accuracy at low energies, i.e., with 
increasing Af. 

In principle, QMC, ED, and NRG are exact methods, but they require an extrap- 
olation, i.e., the discretization of the imaginary time At — » (QMC), the number 
of lattice sites of the respective impurity model n s — > oo (ED), or the parameter 
for logarithmic discretization of the conducting band A — > 1 (NRG), respectively. 

In the present paper, we will not pre sent further, details of t hese, methods and 
refer the reader to the literature: IPT@>B0, NCA@, and QMC00. In the con- 
text of LDA+DMFT we refer to the computational schemes to solve ±he DMFT 
equations discussed above as LDA+DMFT(X) where X=IPTQ, NCA@, QMC0 
have been investigated in the case of Lai—pSr^TiOs. The same strategy was for- 
mulated by Lichtenstein and Katsnelsonlla as one of their LDA++ approaches. 
Lichtenstein and Katsnelson applied LDA+DMFT (IPT)E2, and were the first to 
use LDA+DMFT(QMC)E3, to investigate the spectral properties of iron. Liebsch 
and Lichtenstein also^x>plicd LDA+DMFT(QMC) to calculate the photoemission 
spectrum of S^RuOjS 

2.5. Self-consistent LDA+DMFT 

In general, the DMFT solution will result in a change of the occupation of the 
different bands involved. This changes the electron density p(r) and, thus, results 
in a new LDA-Hamiltonian -Hlda ( |ll| ) since -Hlda depends on p(r). At the same 
time also the Coulomb interaction U changes and needs to be determined by a new 
constrained LDA calculation. In a self-consistent LDA+DMFT scheme, i?LDA and 
U define a new Hamiltonian ( |l8|) which again needs to be solved within DMFT, 
etc., until convergence is reached: 

DMFT 

p(r) H LDA , U n ilm -p(r) 

| | (22) 



Without Coulomb interaction (U = 0) this scheme reduces to the self-consistent 
solution of the Kohn-Sham equations. 

2.6. Simplifications for transition-metal oxides 

Many transition metal oxides are cubic perovskites, with only a slight distortion of 
the cubic crystal structure. In these systems the transition metal d-orbitals lead to 
strong Coulomb interactions between the electrons. The cubic crystal-field of the 
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oxygen causes the d-orbitals to split into three degenerate t^g- and two degenerate 
e g -orbitals. This splitting is often so strong that the tig- or e 9 -bands at the Fermi 
energy are rather well separated from all other bands. In this situation the low- 
energy physics is well described by taking only the degenerate bands at the Fermi 
energy into account. Without symmetry breaking, the Green function and the self- 
energy of these bands remain degenerate, i.e., G q i m , q >i> m '{z) = G(z)S q i„ ltq >i> m > and 
£>qlm,q'l'm>(z) = E(z)6qi m ,q>i> m ' f° r I = Id and q — q d (where l d and q d denote the 
electrons in the interacting band at the Fermi energy). Downfolding to a basis 
with these degenerate g^-Z^-bands results in an effective Hamiltonian fl^DA ( wnere 
indices I — Id and q — qd are suppressed) 

GwM = r^Jd 3 k[(w-'S(u))5 m , m >-(H^(y)) n , mf ]- 1 . (23) 

Due to the diagonal structure of the self-energy the degenerate interacting Green 
function can be expressed via the non-interacting Green function G°(u>): 

G(u>) = G°(u>-SH)= [de N ° (e) . (24) 

J u-Y,(u)-e 

Thus, it is possible to use the Hilbert transformation of the unperturbed LDA- 
calculated density of states (DOS) N°(e), i.e., Eq. (p|), instead of Eq. @. This 
simplifies the calculations considerably. With Eq. (p4J ) also some conceptional sim- 
plifications arise: (i) the substraction of -Hlda m ( M ) on ly results in an (unimpor- 
tant) shift of the chemical potential and, thus, the exact form of -Hlda ^ s irrelevant; 
(ii) Luttinger's theorem of Fermi pinning holds, i.e., the interacting DOS at the 
Fermi energy is fixed at the value of the non-interacting DOS at T — within a 
Fermi liquid; (iii) as the number of electrons within the different bands is fixed, the 
LDA+DMFT approach is automatically self-consistent. 

In this context it should be noted that the approximation Eq. fl24| ) is justified 
only if the overlap between the ti g orbitals and the other orbitals is rather weak. 

3. An example: Lai_ a; Sr :E Ti03 

The stoichiometric compound LaTi03 isiixubic perovskite with a small orthorhomz 
bic distortion (Z TA_— O — Ti w 155°)l23 and is an antiferromagnetic insulatodld 
below Tjv = 125 Kcil. Above Tjv, or at low Sr-doping x, and neglecting the small 
orthorhombic distortion (i.e. considering a cubic structure with the same volume), 
LaTi03 is a strongly correlated, but otherwise simple paramagnet with only one 
3d-electron on the trivalent Ti sites. This makes the system a perfect trial candidate 
for the LDA+DMFT approach. 

The LDA band-structure calculation for undoped (cubic) LaTi03 yields the DOS 
shown in Fig. |l| which is typical for early transition metals. The oxygen bands, 
ranging from —8.2 eV to —4.0 eV, are filled such that Ti is three-valent. Due to the 
crystal-field splitting, the Ti 3d-bands separates into two empty e g -bands and three 
degenerate i2g-bands. Since the i2g-bands at the Fermi energy are well separated 
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Energy, eV 

Figure 1: Densities of states of LaTiC>3 calculated with LDA-LMTO. Upper fig- 
ure: total DOS; lower—figure: partial t^g (solid lines) and e g (dashed lines) DOS 
[reproduced from RefE3]. 

also from the other bands we employ the approximation introduced in section 2.5 
which allows us to work with the LDA DOS [Eq. instead of the full one-particle 
Hamiltonian [Eq. @]. In the LDA+DMFT calculation, Sr-doping x is taken into 
account by adjusting the chemical potential to yield n = 1 — x = 0.94 electrons 
within the £2g-bands. There is some uncertainty in the LDA-cakuilated Coulomb 
interaction parameter (7^4 — 5 eV (for a discussion see Ref£3) which is here 
assumed to be spin- and orbital-independent. In Fig. |[ results for the spectrum 
of La .94Sro.o6Ti0 3 as calculated by LDA+DMFT(IPT, KCA, QMCJ for the same 
LDA DOS at I « 1000 K^and U = 4 eV are compared!!. In RefH the formerly 
presented IPTlI and NCAE3 spectra were recalculated to allow for a comparison 
at exactly the same parameters. All three methods yield the typical features of 
strongly correlated metallic paramagnets: a lower Hubbard band, a quasi-particle 
peak (note that IPT produces a quasi-particle peak only below about 250K which 
is therefore not seen here), and an upper Hubbard band. By contrast, within LDA 
the correlation-induced Hubbard bands are missing and only a broad central quasi- 
particle band (actually a one-particle peak) is obtained (Fig. |l|). 

While the results of the three evaluation techniques of the DMFT equations 
(the approximations IPT, NCA and the numerically exact method QMC) agree 
on a qualitative level, Fig. || reveals considerable quantitative differences. In par- 
ticular, the IPT quasi-particle peak found at low temperatures (see right inset of 
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Figure 2: Spectrum of Lao.94Sr .o6Ti0 3 as calculated by LDA+DMFT(X) at T = 
0.1 eV (w 1000 K) and U = 4 eV employing the approximations X=IPT, NCA, and 
numerically exact QMC. Inset left: Behavior at the Fermi level including the LDA 
DOS. Inset right: X=IPT and NCA spectra at T = 80 K [reproduced from RefSj. 

Fig. |||) is too narrow such that it disappears already at about 250 K and is, thus, 
not present at T « 1000 K. A similarly narrow IPT quasi-particle peak was found 
in a three-band model study with Bethe-DOS by Kajueter and KotliaiS. Besides 
underestimating the Kondo temperature, IPT also produces notable deviations in 
the shape of the upper Hubbard band. Although NCA comes off much better than 
IPT it still underestimates the width of the quasiparticle peak by a factor of two. 
Furthermore, the position of the quasi-particle peak is too close to the lower Hub- 
bard band. In the left inset of Fig. [| the spectra at the Fermi level are shown. At 
the Fermi level, where at sufficiently low temperatures the interacting DOS should 
be pinned at the non-interacting value, the NCA yields a spectral function which 
is almost by a factor two too small. The shortcomings of the NCA-results appear 
to result from the well-known problems which this approximation scheme encoun- 
ters already in the single-impurity Anderson model at low temperatures and/or 
low frequencieseJei Similarly, the deficiencies of the IPT-results are not entirely 
surprising in view of the semi-phenomenological nature of this approximation, es- 
pecially for a system off half filling. This comparison shows that the choice of the 
method used to solve the DMFT equations is indeed important, and that, at least 
for the present system, the approximations IPT and NCA differ quantitatively from 
the numerically exact QMC. 

Photoemission spectra provide a direct experimental tool to study the elec- 
tronic structure and spectral properties of electronically correlated materials. A 
comparison of LDA+DMFT(QMC) at 1000 K0 with the experimental photoemis- 
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Energy (eV) 

Figure 3: Comparison of the experimental photoemission spectrumEl, the LDA 
result, and the LDA+DMFT(QMC) calculation for La .94Sr .06TiO 3 (i.e., 6% hole 
doping) and different Coulomb interaction U — 3.2, 4.25, and 5 eV [reproduced 
fromRefB]. 



sion spcctrunO of Lao.94Sro.o6Ti03 is presented in Fig |[ To take into account the 
uncertainty in £/§!, we present results for U = 3.2, 4.25 and 5 eV. All spectra are 
multiplied with the Fermi step function and are Gauss-broadened with a broad- 
ening parameter of 0.3 eV to simulate the experimental resolution^]. LDA band 
structure calculations, the results of which are also presented in Fig. ||, clearly fail 
to reproduce the broad band observed in the experiment at 1-2 eV below the Fermi 
energyH. Taking the correlations between the electrons into account, this lower 
band is easily identified as the lower Hubbard band whose spectral weight origi- 
nates from the quasi-particle band at the Fermi energy and which increases with 
U . The best agreement with experiment concerning the relative intensities of the 
Hubbard band and the quasi-particle peak and, also, the position of the Hubbard 
band is found for U = 5 eV. The value U = 5 <e.V is still compatible with the ab ini- 
tio calculation of this parameter within LDAt3. One should also bear in mind that 
photoemission experiments are sensitive to surface properties. Due to the reduced 
coordination number at the surface the bandwidth is likely to be smaller, and the 
Coulomb interaction less screened, i.e., larger. Both effects make the system more 
correlated and, thus, might also explain why better agreement is found for U = 5 
eV. Resides that, also the polycrystalline nature of the sample, as well as spin and 
orbitalo fluctuation not taken into account in the LDA+DMFT approach, will lead 
to a further reduction of the quasi-particle weight. 
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4. Conclusion 

In this paper we described the set-up and presented results of the computa- 
tional scheme LDA+DMFT which merges two non-perturbative, complementary 
investigation techniques for many-particle systems in condensed matter physics. 
This approach allows one to perform ab initio calculations of strongly correlated 
electronic systems. Using the band structure results calculated within local density 
approximation (LDA) as input, the missing electronic correlations are introduced 
by dynamical mean- field theory (DMFT). On a technical level this requires the 
solution of an effective self-consistent, multi-band Anderson impurity problem by 
some numerical method (e.g. IPT, NCA, QMC). Comparison of the photoemis- 
sion spectrum of Lai-^Sr^TiOs calculated by LDA+DMFT using IPT, NCA and 
QMC reveal that the choice of the evaluation method is of considerable importance. 
Indeed, only with the numerically exact QMC quantitatively reliable results are ob- 
tained. The results of the LDA+DMFT(QMC) approach are found to be in very 
good agreement with the experimental photoemission spectrum of Lao.94Sr .o6Ti03. 

The LDA+DMFT(QMC) scheme will provide a powerful tool for all future ab 
initio investigations of real materials with strong electronic correlations. In particu- 
lar, LDA+DMFT is the only existing ab initio approach which is able to investigate 
correlated electronic systems close to a Mott transition, as well as heavy fermion 
and /-electron materials. The physical properties of such systems are character- 
ized by the correlation-induced generation of small, Kondo-like energy scales which 
require the application of genuine many-body techniques. 

With LDA+DMFT the band structure and model Hamiltonian communities 
which essentially lived separate lives so far are finally able to join forces. In the 
future both communities will need each other's input. Indeed, without DMFT the 
LDA will be limited to the investigation of weakly correlated systems, and without 
the LDA-calculated band structure the many-body approach will be restricted to 
the study of more or less oversimplified models. 
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